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Abstract. In the field of radial basis functions mathematicians have 
been endeavouring to find infinitely differentiable and compactly supported 
radial functions. This kind of functions are extremely important for some 
reasons. First, its computational properties will be very good since it's com- 
pactly supported. Second, its error bound will converge very fast since it's 
infinitely differentiable. However there is hitherto no such functions which 
can be expressed in a simple form. This is a famous question. The purpose 
of this paper is to answer this question. 
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1 Introduction 



In this paper we raise a new function 



— a 



$(x) : 



gTHNP if || x || < 1 

if \\x\\ > 1 



(1) 




where a > is a constant and denotes the Euclidean norm of x e i? n . 
The infinite differentiability of $ can be shown by routine investigation, and 
we omit it. What's difficult is to show its positive definiteness. In order 
to make it useful in scattered data interpolation, this problem has to be 
overcome. This is the main theme of the next section. 

2 Positive Definiteness 

It's well known that positive definiteness is just conditional positive definite- 
ness of order zero. Here we adopt the definition of [1] for conditional positive 
definiteness. Before showing that the map defined in (1) is positive definite, 
let's define 



where a > is a constant. In other words, (f>(\\x\\) = §(x) for x <E R n . In 
order to make this paper easier to understand, let's analyze the derivatives 



• f t (f>(Vi) = e- au (-a)u 2 

• &<f>(Vi) = e~ au [a 2 u 4 - 2cm 3 ] 

• ^4>{V~t) = e- au [-a 3 u 6 + Qa 2 u 5 - 6cm 4 ] 

• £<j>(y/i) = e- au \a A u 8 - 12a V + 36a 2 u 6 - 24cm 5 ] 

• ^4>{Vt) = e- au [-a 5 u 10 + 20cAr> - 120a 3 u 8 + 240aV - 120cm 6 ] 

• ^(p(Vt) = e- au [a 6 u 12 - 30a 5 u u + 300a 4 w 10 - 1200aV + 1800a 2 u 8 - 





of (t>(\/i) of order up to eight. 
It's easily seen that with u 



= j-^, we have ^ = u 2 for u > 1 and 



720cra 7 ] 
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• = e -^_ a 7 u u + 42a 6 u 13 - 630a V 2 + 4200a V 1 - 12600a V° + 
15120aV - 5040cm 8 ] 

• £t<j>(Vi) = e- QU [a 8 u 16 -56aV 5 +1176aV 4 -11760aV 3 +58800a 4 w 12 - 
141120a 3 « n + 141120aV° - 40320aw 9 ]. 

By induction, for j > 1, 

(P 

—<j>(Vi) = e- au [(-l) j a j u 2j + (-l) j - 1 a j - 1 (x>- 1 u 2j - 1 + ■■■ + (-l) ai au j+1 } 

(XV 

where a±, • ■ • , aj„i are nonnegative integers. Obviously 

(-l) J ^(v^) > 

for all j > 1 and n > l(i.e. < t < 1) as long as a is large enough. In fact 
we have the following more accurate lemma. 

Lemma 2.1 Assume that «i = and ojj+i = a, + 2j /or j = 1,2,3, • • •. 
Then for any j > 1, 

<-i4^)>o 

for all t > if a > ctj . 

Proof. Note that if ctj is defined as in the lemma, then aj is just the exponent 
of u in the first term of f tj _ l (f)(Vi) plus the absolute value of the integer 
coefficient in the second term of <fi(V~t)- In turn the absolute value of the 
integer coefficient in the second term of 4rj=r(t>(Vi) is the exponent of u in 
the first term of J^ 7 0(y / t) plus the absolute value of the integer coefficient 
in the second term of -^r^&iVi)- Our lemma thus follows immediately by 
simple induction. 

Theorem 2.2 (a)Suppose a > 12. Then $(a;) in (1) is positive definite on 
R d for 1 < d < 4. (b)Suppose n > 4 and a > aj for j = l 1 ^] + 4, where 
[A] denotes the least integer not less than X. Then $(a;) in (1) is positive 
definite on R d for 1 < d < n. 

Proof. This is an immediate result of Theoreml of [1], Theorem2.2 of [2] and 
Lemma2.1 of this paper. 
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Remark: In this theorem we treat R d for 1 < d < 4 independently to 
make it easy for use. In fact, if we require $ be positive definite on R d 
for 1 < d < 2 only, then a > 6 is sufficient. The first values of aj are 
0, 2, 6, 12, 20, 30, 42, 56, 72, 90, • • • ; for j = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, • • •, respec- 
tively. It's easily seen from our deduction of -^j(p(y/i) that 

(-l)'^(Vf) > 

for j — 1, 2, • • • , 8 if a > aj, and t > 0. Thus we already checked that 
is positive definite on R d for 1 < d < 12 if a > 56. For higher dimensions, 
the values of the necessary aj can be calculated recursively. It's a bit more 
tedious, but very easy. 
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